Asymptotically exact mean field theory for the Anderson model including double 

occupancy 



J. Holm, R. Kree, and K. Schonhammer 
Institut fur Theoretische Physik, Universitdt Gottingen, 
Bunsenstrafle 9, D-3400 Gottingen, Germany 
(March 1993) 

The Anderson impurity model for finite values of the Coulomb repulsion U is studied using a 
slave boson representation for the empty and doubly occupied /-level. In order to avoid well known 
problems with a naive mean field theory for the boson fields, we use the coherent state path integral 
representation to first integrate out the double occupancy slave bosons. The resulting effective 
■ action is linearized using two-time auxiliary fields. After integration over the fermionic degrees of 

freedom one obtains an effective action suitable for a l/A^-expansion. Concerning the constraint 
ON . the same problem remains as in the infinite U case. For T — > and Nf — > oo exact results for the 

ground state properties are recovered in the saddle point approximation. Numerical solutions of the 
saddle point equations show that even in the spindegenerate case Nf = 2 the results are quite good. 
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I. INTRODUCTION 



The Anderson impurity model (AIM)0 was proposed 
' to gain insight into the physics of dilute magnetic alloys. 
^ . Variations of this model describe e.g. thermodynamic, 
transport and spectroscopic properties of Ce mixed va- 
lence compounds surprisingly wellocl. The strong on-site 
/-repulsion U makes these compounds typical examples 
of strongly correlated electron systems. 

Various reviews exist on the theoretical approaches to 
| obtain exact or approximate solutions for thepjhermody- 
namics and correlation functions of the AIMIJB. While 
exact solutions for thermodynamic properties are known 
for special cases using the Bethe-AnsatzEl, the progress 
towards controlled approximations for spectral functions 
was mainly due to large AU expansions, where Nf is the 
5 ■ degeneracy of the f-leveJ3'Q 



degeneracy of the /-level 

The two special cases for which the exact thermody- 
• • i namics can be obtained via the Bethe-Ansatz are the 
infinite U limit (U — oo,Nf arbitrary ) and the spin- 
degenerate symmetric case (2e/ + U — 2[i,Nf = 2) 
' where n is the chemical potential which will be chosen 
as the zero of the one electron energies in the follow- 
ing and £f is the bare /-level energy. The one-particle 
spectral functions which determine photoemission and in- 
verse photoemission spectra are rather well understood 
in the infinite U limit for arbitrary temperatures and de- 
generacies down to Nf = 2. For temperatures of the 
order of the Kondo temperature Tk or smaller and /- 
occupancies close to one (spinfluctuation limit) there is a 
sharp " Kondo- resonance" slightly above the Fermi en- 
ergy and a broad resonance near Sf. The knowledge 
about the position, the width and the weight of the 
Kondo resonance comes from three different approaches. 
It was first obtained for zero temperature using a varia- 
tional ansatz for the grourujLstate and a l/Nf classifica- 
tion of intermediate statestrQ. The noncrossing approx- 
imation (NCA) which can be derived by different t 
niques works well for all temperatures except T<T J 
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In the third approach a functional integral (FI) for the 
partition function is expressed in terms of an effective ac- 
tion for the slave boson, which is introduced to handle the 
constraint of at most single occupancy of the /-leve£j~£3. 
If the FI is evaluated in the saddle point approximation, 
corresponding to a simple mean field approximaliguJbr 
the slave boson, the exact ground state energyB0i29 is 
recovered in the limit Nf — > oo. Gaussian fluctuations 
around the saddle point yield (1/Nf)— corrections in the 
very low temperature Fermi liquid regimeEj. Due to the 
occupancy constraint the saddle point expansion is not a 
(1/Nf)— expansion in the strict sense. This problem will 
be addressed later in this paper for finite values of U. 

In the spindegenerate symmetric case the Kondo peak 
is at the FermiJej/el due to the particle-hole symmetry 
of the problema'B. For Nf > 2 this symmetry is no 
longer given even for symmetric parameters 2e f + U = 
and a symmetric band, and the position of the Kondo 
peak haSptipt-been reliably calculated except in the large 
Nf Iimitil3il3. Also the position of the Kondo peak 
for nonsymmetric parameters and Nf = 2 is not yet 
well understood. For not too low Kondo temperatures 
the only reliable results are from quantum Monte Carlo 
calculationstL3. The spectral function for the Nf = 2- AIM 
has recently gained renewed attention as its calculation 
is required as an intermediate step in the ex-ajt-solution 
of the Hubbard model in infinite dimensionst-J~li3. 

For the finite U AIM and arbitrary degeneracy Nf 
there exists a variational approacbllij at T = and gen- 
eralizations of the NCAoliS pjwhich go over to the cor- 
responding U — oo theoriesujafl in the large U limit. A 
naive mean field theory for the slave bosons for the empty 
/-level as well as the doubly occupied /-level does not re- 
cover the exact ground state energy in the limit Nf — > oo, 
as it completely misses to obtain the mathematical struc- 
ture of an integral equation-.f£pJinite C/EJ in contrast 
to an transcendental equationLrQEl for U — oo. Similar 
deficiencies of mean field theories for double.-ajCcupancy 
slave bosons occur for the Hubbard modelE^cS. 

It is the aim of this paper to present a generalized 
mean field theory for the finite U Anderson model 
in the Z ,/ 1 and / 2 subspace for arbitrary degeneracy 
Nf, which leads to the exact ground state energy in 
the limit Nf — > oo. Starting from a coherent state 
path integral we integrate out the double occupancy 
bosons, which leads to four - fermion terms. Using a 
Stratonovich-Hubbard transformation which introduces 
two-time auxiliary fields an effective action appropri- 
ate for a large Nf expansion is derived, apart from the 
problem with the occupancy constraint, similar to the 
infinite U caseOJ. We study stationary solutions of the 
saddle point equations for which the additional auxiliary 
fields depend on the time difference only. This leads 
after Fourier expansion to a system of coupled nonlinear 
equations for the Fourier components. Analytic contin- 
uation from the Matsubara frequencies to the real axis 
leads to a system of integral equations, which for finite 
Nf has to be solved numerically. In the limit Nf — > oo 
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the equations simplify considerably and for T — > we 
exactly recover the integral equation which determines 
the ground state energy^. 

In section II we present the model and introduce the 
coherent state path integral formulation. An exact ef- 
fective action for the empty /-level slave boson and the 
two-time auxiliary fields is derived. In section III the 
corresponding saddle point equations are derived. The 
simplified equations in the limit Nf — > oo are discussed. 
The numerical solutions for finite Nf are presented in 
section IV. Some results of ref.EJ and details of various 
calculations are described in appendices. 



II. THE MODEL 



We cons, 
Hamiltoniar' 1 ' 



the finite U Anderson impurity 



T I 1 J j V.tl 



v=l 

= H + V + U 

where we have introduced a combined index, is, for the 
orbital and spin degeneracies. Spin-orbit and multiplet 
splittings are neglected. The first term describes the con- 
duction states, with energy e, and the second term the 
/-level with the (bare) energy £/. The third term leads 
to a hopping between these states where the V(e) can 
be chosen real and the last term describes the Coulomb 
interaction between the /-electrons. The Hamiltonian 
is studied in the subspace of at most double occupancy 
of the /-level i.e. in the f ,/ 1 and f 2 manifold. This 
restriction is incorporated using slave bosons. 

In the infinite U model Colemana has introduced a 
single slave boson as a counting device for the number of 
/-electrons. For the finite U model there are Nf(Nf — 
l)/2 different doubly occupied /-states. Therefore, we 
introduce N f (N f - l)/2 additional "heavy bosons"liI If 
| {0} > denotes the ground state of the "light" and heavy 
bosons we make the correspondence 



}>®|{0}> (2) 



A* > v : /+/+ | K} >- 

where | {n E } > denotes an arbitrary band state and 
hz v is the creation operator of a heavy boson. To remove 
the restriction fj, > v it is useful to define boson operators 
for n < v by h+ v := and = 0. 

We now construct a Hamiltonian acting on the product 
space of electrons and bosons which is equivalent to H in 
the {Z ,/ 1 ,/ 2 } subspace 



(1) 



H = H + Y / (^f + U)h+„h^ 



de 



N f 

E 



N f 



V(e)f+c ev b+ V(e)h+,f vCcM +h.c. 



(3) 



Here Ho is formally the same as in Eqn. (|l|). The 
Hamiltonian H commutes with the " charge operator" Q 
defined as 
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Q = b+b+Y,rtu + Y, h ^ h ^ ( 4 ) 



Therefore the subspaces Fq with definite integer value 
Q can be treated separately. In the subspace F\ the 
Hamiltonian H is equivalent to the Anderson Hamilto- 
nian Ha in the subspace of J ,/ 1 and f 2 states. The 
partition function is therefore given by 



Z A = Tne-? 11 = Z 1 



(5) 



where the trace is restricted to the subspace of states 
with Q = 1. This projection to the Q = 1 subspace can 
be expressed as a contour integral in a complex A-plane 
(Z = Z A ) 



Z = A [ d\e' 3X Tr(e-^ x A 
2ni J c \ J 



(6) 



where H{\) = H + XQ and the contour runs from 
\r — in/ f3 to Xfj + in/13, where the real \r is arbitrary. 
Since the trace has now to be performed over the full 
Hilbert space of fermions and bosons Z can be expressed 
as a coherent state functional integral 



(3d\ px 



z = 


L 






where 






Nj 


L (r) = 


E 



- — e pA / DfDcDbDhexp- 
2ni 



(L (t) + ii(r))dr 



+b*(r) 



dec ev (T) 



d ^ 

0~r +X 



d_ 



d_ 



■e f + X 



&w + E^( 7 



— + 2e f + U + \ 

OT 



Mr) 



and 
Li(t) 



]T / ds[V(s)f v (r)c El/ (T)b(T)+V(e)b*(r)c ev (r)f v (r)] 
- J2 / * WW/.WCcW + V(s)^(t)/»(t)VW] 



with r-dependent Grassmann fields c, c and /, / and 
complex fields b and h^ v . 

As the Grassmann fields enter the action S = 
Jq c?tL(t) with L(t) = L (t) + L\{t) quadratically it 
would be possible at this point to eliminate all fermionic 
variables by tracing over in the usual way. As we want to 
construct a mean field theory which becomes asymptoti- 
cally exact in the large degeneracy and small temperature 
limit this procedure turns out not to be successful. The 
saddle point equations for the resulting effectivje-.action 
S e f f (b, h) do not produce the exact T — results^ in the 
large degeneracy limit. This failure can be related to the 
fact that the term Lih(r) in L\(t) is not well suited for a 
large Nf saddle point expansion,because it consists of Nf 
terms in which different orbital indices /i, v are coupled. 
To bring this term in the action into a suitable form, 



we proceed differently, in a way which in the first step 
produces a seemingly more complicated effective action 
involving " four-fermion" terms. Factorizing these terms 
with a Stratonovich-Hubbard transformation leads to an 
interaction term in the effective action of the form needed 
for a large Nf expansion. 

In the first step the heavy bosons are integrated 
out. The Gaussian integrals are easily performed and 
yield using h vli = —h^v in the interaction term Li^(r) 





r r fi ' 


j Dh exp < 






[ Jo 



E K»( T )(]£ +2s f + U + A) fe F (r) + L lh (r) 

4(A) exp { - V / dede' / dr \ dr'V{e) (c^(t)U(t) - c^(r)/ Al (r)) 
„^,jJ Jo Jo 



G° h (r - t>) (U(r')c^(r') - f,{r')c^{r')) V{e')}, (10) 



where 

4(A) = (1 - exp [-0(2e f + U + x)]y N f( N f- 1 )/ 2 is the 
unperturbed partition function for the heavy bosons and 

G° h (r - r') = -(r | {d T + 2e f + U + A)" 1 | r') (11) 

is the corresponding unperturbed propagator with the 
usual bosonic boundary conditions. 

The exponent on the rhs of Eqn. ([To]) can be written 
in a form suitable for a Stratonovich-Hubbard transfor- 
mation introducing 



Vc^{t) = J deV(s)c^(T) (12) 

with V = (J* | V(e) | 2 de) 1/2 and using 

£Mt)/„(t) " 5u(r)U(r))(fu(r')c,{r') - f,{r')c v {r')) (13) 

(E ^(r) C(I (r'))£ /„(r)/„(r')) + £ ^(r)/,^))^ /„(r)c„(r')). 



As both terms on the rhs have a simple product form 
we can use the identity 

e CD _ A J dfiea dlma c _ A]a] 2_ VJ(aC+a , D) ^ 

valid for commuting variables C and D and i?e(A) > 
0. The auxiliary fields which have to be introduced to 
linearize the four-Grassmann terms depend on two time 
variables 



exp j -V 2 (J2 ^(tMt'))(]T Ut)U{t'))GI{t - t')\ 

L - J dn(F(r, r')) exp j ^ - F(r, /) V £ ^(r)c^(/) - F* (r, t> £ / M (r)/ M (r') 



where d/i(F) = dReFdlmF/n. The second term in 
Eqn. (ph yields 
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cxp (-^ 2 (]T c m (t)/ m (t'))(£ Ut)c u (t'))GI{t - t>)\ 

L fj, v ) 



Gib 



d M (X(r,r'))exp|^^ - V £(X(r, t'^t)/^') + X*(r, r')/,(r)c M (/)) 



The linearization produces additional terms bilinear in 
the Grassmann variables. The first Grassmann term in 
the exponent on the rhs of Eqn. ( |l5| ) corresponds to a 
separable perturbation ~ | c, fi >< c, fj, | , where 



c,(i>:= i y cfeV(e) | e,ju > 



(15) 



are linear combinations of conduction states localized 
around the impurity. The other term in the exponent on 
the rhs of Eqn. (|l5|) leads to a fluctuating /-level position 
Ef+\— > £/+A+V^F*(r, r'). The Grassmann terms in the 
exponent on the rhs of Eqn. (2.14) modify the hopping 
strength V&(t)<J(t - r') -> F(&(t)£(t - r')) + X*(t, t')). 

At the next step the Grassmann variables are formally 
integrated out 

Z = J ^ e 0X Z$(\) J DbDFDXe~ Scff (16) 

where the effective action S e f / is given by 

S eff= I b {T)(d T + X)b{T)dT - I dn I dr 2 — ^ 



(i 



o 







iV/Tr In [9(E) l e + l T ®h +hi] (17) 



and the trace runs over the space of antiperiodic func- 
tions on the interval [0, 0\ as well as the space of single 
electron states for a single channel. The operators ho and 
hi are given by 

h = j e\e><s\de + {e f + X)\fxf\ (18) 

hi = V(F* | / >< / | +F | c >< c | (19) 
+ (&*+*) | cx/| +(& + **) | /><c|) 

The label of the electronic channel has been supressed 
as each channel gives the same contribution leading to 
the factor Nf in front of the trace. Scaling all Bose fields 
by a factor y/Nf i.e. b — ► */ 'NfB etc. the effective action 
becomes proportional to Nf and a well defined large Nf 
expansion could be carried out, if, not the same problem 
known from the infinite U casetil due to the constraint 
integral over A would show up. The exponent [3X is not 
proportional to Nf which leads.to a similar behaviour as 
discussed in detail for U = ooEil. 

Nevertheless we treat the functional integral in Eqn. 
( |l6| ) by a saddle point approximation (without scaling of 
the Bosc fields). 
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III. THE SADDLE POINT EQUATIONS 



To derive the saddle point equations (SPE) the 
Bose fields and their complex conjugates (F, F*etc.) are 
treated as independent variables and we replace F* — > F 
etc. to indicate that the saddle point values of F and F 
etc. are in general not complex conjugate fields. 

The functional derivatives of the effective action S e f / 
involve single electron propagators which we abbreviate 
as = f or c) 



G 13 {t,t') = -(t |< i | [d ® l e + l T <8> h + h^- 1 \j>\r') (20) 



which obey the usual fermionic boundary conditions. 
As the factor e l3X Z^ (A) is independent of the Bose fields 
the SPE from the derivatives with respect to the Bose 
fields read 

= ^ = (Or + A)6(r) + N f VG fc (r, r) (21) 
= ^ = {-dr + X)b{r) + N f VG cf (r, r) (22) 
»-^^-- * {Ty) +N f VG fc (r',r) (23) 



5X(t,t>) G h (r-T>) 

5S e ff _ X(t. 



■ T 



6X(t,t>) G° h (r-r') 



N f VG cf {r',r) (24) 



«=^ = -^ + »»V> (26, 

The generalization of the static approximation in the 
U = oo case amounts to search for solutions with " time" - 
translational invariance, i.e. the two-time Bose fields de- 
pend on the "time" -difference only, while b and b are 
constant. This yields the equations 





- Xb 


= N f VG fc (0) 




(27) 




-Xb 


= NfVG cf {0) 




(28) 


X(r 


-r') 


= N f VG° h (T- 


r')G fc (r'-r) 


(29) 


X(r 


-r') 


= N f VG° h (r- 


r')G cf {r'-T) 


(30) 


F(r 


-r') 


= N f VG° h (r- 


t')G cc {t> - t) 


(31) 


F(r 


-r') 


= N f VG° h (r- 


r')G ff (r'-r) 


(32) 



In the SPE from the constraint variable A we have also 
to take into account the A-dependence of the integrand 
in dl5| ) not contained in S e ff. As we are interested in the 
low temperature regime the partion function Zq(X) can 
be approximated by 1, but the factor e ,3A has to be kept. 
This yields 



0= ^//-^) (33) 
oX 

b(Ti)b(ri)dr - N f J G ff (r,T + 0)dr - (3 

J (X(T U T 2 )X(T U T 2 ) + {F(n,T2)F(T 1 ,T2))-^(l/G h {T 1 -T 2 ))dT 1 dT2 
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Using Eqn. (§9|-||) this SPE reads ((n f ) ps = 
NfG ff (-0)) 

\b\ 2 +(n f ) ps +A = l (34) 
where A is given by (V 2 = NfV 2 ) 

A = V 2 J d T[N f G fc (-T)G cf (-T)+N f G ff (-T)GU-T)}-^G° h (T) (35) 

As will be shown in appendix C, A presents the prob- 
ability for the double occupancy of the /-level. To solve 
the SPE ( p7[ - p^ ) and (]34|) we write the propagators as 
Fourier series 

G °h ( T - r ') = \ E e- iVm{T - T,) Gl (iv m ) (36) 

m 

GiAr' ~ r) = - e-^'-^G^iuJn) (37) 

h 

where the Matsubara frequencies u m — 2?:m/(3h,m G 
Z are of the Bose type, while the uj n = 2ir(n + 1)/ (3h, n e 
Z are fermionic frequencies. Then Eqns. ( p7| - |3^ ) lead to 
equations for the Fourier coefficients, e.g. 

X(uu n ) = N f V& G° h (iu m )G fc (w m - (38) 

m 

and corresponding equations for the other Bose fields. 
To calculate the Gij(iu n ) one has to calculate the resol- 
vent of the one electron operator h(iuj n ) = Hq + hi(iuj n ) 

Gij(iuJ n ) =< i | (iuJ n - /i(«w n )) -1 | j > (39) 

The inversion is straightforward using a partioning 
technique with the projector P =| / >< / |. With 

G»( Z )e< c |( z -/ 1o m )- 1 |c> (40) 

and 

G° cc (icj n ) = G£(*w»)/(1 - VF(iu n )G™(iu n )) (41) 

the result for the f- propagator is given by (if = Ef + A) 

Gff(icj„) = [ito n -e f - VF(iiu n ) ~V 2 (b + X(iuj n ))G cc (icj n ){b + X(iu n ))}- 1 (42) 

The other resolvent matrix elements are 

G cf (iij n ) = G° cc {tLu n )V(b + X(iu n ))G ff {iLu n ), (43) 
G fc (iuj n ) = G ff (iw n )V{b + X(iu; n ))G cc (iLu n ), (44) 
G cc (iu n ) = G° cc {iiu n ) + G cf (ito n )V(b + X{nu n )) (45) 

Alternatively G cc can be written in the form 

G cc {iuj n ) = (iuj n - if - V F(iu n ))G f f(iuj n )G cc (iu n ) 

(46) 

In the following we study the solution of the SPE with 
b, b different from zero and define 
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x{iuj n ) = X{iuj n )/b (47) 
x(iuj n ) = X(iuj n )/b (48) 



Then the SPE (|27fgg) take the form 



- A = ^ 2 4E G //( iW ™) G °c(^n)(l +X(kj n )) (49) 
n 

x(iuj n ) = y 2 i Gl(iis m )Gff(iv m - za; n )G|? c (zz/ m - (1 + x(iu m - iw n )) (50) 

m 

y_F(zcj n ) = y 2 i ^ Gl(iu m )G^ c (iu m - iw n ) (ii/ m - iio n - if - VF(iv m - iw n ))Gjf{iv m - ico n ) (51) 

m 

VF(iu n ) = Vi-YjGliiv^GffiiVrn-iUn) (52) 

m 

The equations following from (E8h and (pOf ) involve 
x(i(jj n ). Comparing with Eqns. ( [49| ) and (pQj leads to 
the identity x(icu n ) = x(iui n ), i.e. we can drop the addi- 
tional equations. Together with (j3j|) the Eqns. 
determine the Bose fields at the saddle point. 

In order to obtain a first understanding of these equa- 
tions we start with a discussion of the limit Nf —* oo. 
It is easy to see that it is consistent to assume that 
x(iu> n ), VF(iu> n ) and VF(iui n ) are of order unity in the 
limit Nf — > oo with V constant. This implies 

Jy' — #■ oo If) 

Gf f (iu) n ) — ► (iu n - if - VF(iu> n ))~ = G ff (iu n ) 

(53) 

and G cc (iuj n ) — ► G Q cc {iu) n ). This simplifies Eqns. (|5 
and (|5l]) for F(iu n ) and F(iuo n ) 

VF(iu> n ) = V 2 W 



P ~ \w m - (2e/ + U + A)] [iu m - iu> n - if ~ VF(iv„ 

(54) 



VF^ n ) = V^^2- -1 



+ C7 + A)] [(G™(iv m - iu n ))-i - VF{iv m - iu n )\ 

(55) 



which now constitute a closed system of equations for 
a given value of A. At this point it is necessary to discuss 
the strategy to solve the SPE for the Fourier components. 
Guided by results of an expansion in powers of 1/U we 
will assume that an analytic continuation F(ioj n ) — > 
F(z) etc. is possible, with F(z) etc. analytic functions 
except at (parts of) the real axis. As F and F 
are of order 1/U the leading order result for VF(iuj n ) is 
obtained by neglecting VF(iv m —iLu n ) in the denominator 
in Eqn. (p4). Then the Matsubara sum can be performed 
in the usual way by contour integration and one obtains 

VF & n ) = V> n -&' + U + » + K e ' + * (56) 
iu) n -Sf - U 

where n-{e) is the Bose function and f(e) is the Fermi 
function. The leading order result for VF(iuj n ) is 
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J iw„ + e - (2e/ + (7 + A) 

(57) 

To this order the analytical continuation is trivial and 
has the properties discussed above. In the low tempera- 
ture limit the Bose functions at ijj = 2ej + U + A can be 
neglected. At T = the function VF^{z) has a branch 
cut from iu to ijj + B. For the solution of Eqns. ( |5^ , 
|55| ) for arbitrary [7 we now assume the spectral repre- 
sentations (i = c, /) 

G^(i. n )= p^-de (58) 



Performing the Matsubara sums in Eqn. (|54|, |55|) yields 
after analytic continuation 

VF(z) = y, [tifiWW +»-&)) ^ 
J z + e-iu 

vf(z) = J m 

Together with Eqn. (g) and G° cc (z)- 1 = G^(z)- 1 - 
VF(z) this constitutes a system of integral equations for 
the spectral weight functions p°f and p° cc . As the finite 
temperature theory will be discussed later for arbitrary 
Nf we restrict ourselves here to the limit T — > 0. At the 
saddle point value of A one expects P°ff{£) to vanish for 
e < Ef as in the U — oo case0. In addition to the Kondo 
peak slightly above the Fermi energy £p = 0, pfl{£) has 
an f 2 contribution above £jj, as can be inferred form 
Eqn. (|57|). Therefore the integrand on the rhs of Eqn. 
( |59| ) vanishes, i.e. VF(z) = at T = in the limit 
N f -> oo. This leads to G° c (z) = G°°0) and 

G° ff (z) = (z- £ f + f (2e, + U + X- z))- 1 (61) 

with r(z) defined in (A. 6). 

For T -> and TV/ -> oo Eqns. @ and @ then 
simplify to 

./-fl P„ [lUJn-Ef+T^u-lUJn)] [l.U n - £) 

= f ' I I 2 ^ V ? 1 , = 

7-B P „ (*^m - EC/) [il/ n - iw„ - E/ + r(e;7 - M^m + iu n )] [}V m ~ lW n - £) 

In order to perform the Matsubara sums by contour 
integration it is necessary to analytically continue the 
Fourier coefficients x(iui n ). A useful guide is again to first 
calculate x(ito n ) to leading order in 1/U, which amounts 
to iterate Eqn. (|6^). This calculation shows that x(ioj n ) 
can be analytically continued and at T — has a branch 
cut from ijj to ijj + B . To evaluate the Matsubara 
sums in ( p^ ) we therefore assume that x{z) has a spectral 
representation with a branch cut starting at a positive 
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energy value. Then for T — > only the last pole term 
on the rhs of (|6^, |6^) contributes to the contour integral 
and one obtains after analytic continuation 



a=/° \ns)m + m) d£ 

J-B e - e f - X + T(2e f + U + \-e) 



' -b (2£ f + U + X-£-z)[e-£f- X + f(2e f + U + X- e)}° 



If z is chosen on the real axis —B< Rez < the second 
equation is an integral equation for x(e) with e G [-B, 0] 
which together with (^54|) determines A. With the identi- 
fication x(e) <-» c(e) and A <-> —AE these are exactly the 
equations which occur in the. description of the ground 
state in the limit Nf — > ooEil This is discussed in ap- 
pendix A. A direct proof for the relation —A = AE valid 
at T = in the limit Nf — > oo is presented in appendix 
B. 

After this discussion of the limit Nf — > oo we return 
to the full SPE (||-|5|) for finite values of Nf. We use 
the same strategy as for the solution of Equs. ( |54|]55| ) 
and assume that the Greens functions Gij(iuj n ) can be 
analytically continued, and as a function of the complex 
variable z are analytic expect at (parts of) the real axis. 
This implies the spectral representation 

Gij(z) = r ^-de (66) 



Using these spectral representations the Matsubara 
sums in Equs. fl49|-|52|) can easily be performed. If we 
also introduce the functions Rij(z), 



we obtain 



VF(z) = -V^n^Gffieu -z)+R f f(i u - z)] 

VF{z) = -V 2 [n-(iu)G cc (iu - z) + R cc (iu - z)\ 
X(z) = ~N f V[n-(eu)G c f(eu - z) + R cf (su - z) 



Therefore F, F and X are also analytic except at the 
real axis and can also be expressed in terms spectral func- 
tions pp(e), pp(e) and pxje)- Together with the relations 
( ^ - |45| ) the equations (|6^-[70|) provide a set of nonlinear 
integral equations for these spectral functions for given 
values of A and b. In order to determine the saddle point 
values of A and b we in addition have to use equation ( p!9| ) 
and the constraint relations ( |34| ) which was irrelevant in 
the limit Nf — > oo. Using the spectral representation 



( pq ) the probability for double occupancy A (35) can be 
expressed in terms of the spectral functions. For temper- 
atures ksT <C §o one obtains 



The numerical solutions of the coupled system of equa- 
tions is described in the next section. 
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IV. NUMERICAL SOLUTION OF THE SADDLE 
POINT EQUATIONS 



In this section we present an outline of our numerical 
procedure to solve the SPE for a given band density of 
states, which determines the unperturbed band propaga- 
tor G™ . The equations ( |68| - |7C| ) provide relations between 
the spectral functions pk{s) with K S {F, F, X} and the 
spectral functions pij (e) with i, j € {/, c}. For fcgT <C iu 
they read 

PF (e) =N f Vpf f (E u -e)f(e u -e) (72) 
p P (e) =N f Vp cc (£u-e)f(eu-e) (73) 
p x (e) = -NfVpfdiu - e)f(iu - e) (74) 

We solve the SPE iterativly using e.g. the results for 
Nr — » oo discussed in section III as the starting values 
for the pk, A and b. From the pk we calculate the func- 
tions K(e + iO), which via equations ( ^2| - ^5| ) determine 
the propagators Gy (e + i O) and their spectral functions 
Pij(e). Using equations ( |72| - |74| ) we obtain new spectral 



functions pif(e). The constraint relation (34) involves 
the pseudo /-occoupancy (n/) ps 

(n f ) ps =N f J Pf f (e)f(e)de (75) 

and the probability for double accupancy A, which we 
evaluate using @ and <^^) 

a /"-j ,7 , Pf{£)Pf{z') + Px{e)px{e') . . , , . 
J dede {eu-e- e') 2 ^ £ >^ £ > ( 76 ) 

In order to obtain a numerically stable method to de- 
termine new values of A and b we treat equations j3^ ) 
and d|) 

1 -(»/)„,- A - ft 2 = (77) 
Xb + N f V f p fe (e)f(e)de = (78) 



as a two dimensional system of equations and solve it 
by Newton's method. We then go to the beginning of 
the loop and terminate the iterative procedure when the 
relative change of all propagators is less then 1CP 4 and 
A and b solve (0) and ([TSh with an error less then 10~ 3 . 
For the functions pk{e) we usually used a 500 points 
logarithmic mesh, which was implemented by logarithmic 
submeshes to resolve characteristic features such as band 
edges. 

In our calculations we used a semielliptical band den- 
sity of states 2VB 2 - e 2 /nB 2 with B = 6. The nu- 
merical results presented in this paper are restricted to 
the zero temperature limit. Compared to the U = oo 
limit the spectral functions pk corresponding to the two 
time auxiliary fields provide a new input for the calcu- 
lation of the psuedo /-Greensfunction, which is given 
by Gff. In figure [l] we show the spectral functions 
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Pf, pp and px for Nf — 14, £/ = —2.5, U = 8 and 
A = 2V 2 /P> = 2.67. The spectral functions are zero for 
energies less then = 2ef + A + E/ w 6.25 which is de- 
termined by Asp ~ 3.25. In the following the subscripts 
SP will be dropped. These spectral functions together 
with A and b determine the pseudo /-spectral function 



(Pff)ps = Pff defined in (|6q ) which is shown in fig |^. 
The main spectral feature is the "Kondo - peak" slightly 
above the chemical potential taken as the zero of energy. 
As can already be seen in the limit Nf — > oo ( |53| ) and 
discussed in detail intj the peak position ist not given by 
Ef + A as (for a flat band density of states) in the U = oo 
limit. There is additional spectral weight in the interval 
[iu,£u + B] which is magnified in the inset of the fig- 
ure. This is f 2 - weight is very small for the pseudo / - 
spectral function, while it is large for the real / - spectral 
function as discussed in referencesEj andtil. The real / - 
spectral function (pff) re al m the- saddle point approxima- 
tion consists of one contribution \b\ 2 (pff(£)) ps , i.e. the 
Kondo peak position in (pff)reai is determined by the 
peak position in (pff) ps - There are additionalspterms in 
the expression for the real / - Greens functionlU which 
are responsible for the large f 2 - weight in {pff) re al- 

The spectral functions p C c(s) and Pfcis) are shown in 
fig|3|. This dip in Pcc(s) at the position of the Kondo peak 
is readily understood from Eqn. (46). The occupancies 
(P n ) for the parameter values used in fig. |l|-||are (P ) = 
0.144, (Pi) = 0.792, (P 2 ) = 0.063. This corresponds to 
the real / - occupancy (nf) rea i = (Pi) + 2(P 2 ) = 0.92. 

The results shown so far correspond to Nf = 14, i.e. 
the large / - degeneracy obtained by neglecting spin - 
orbit and crystal field splitting. As there are no exact re- 
sults available for the model for arbitrary Nf and U the 
deviation of the the SPE results from the exact ones can- 
not be judged. A very stringent test of the quality of the 
saddle point results can be made for the most interesting 
value Nf = 2 of the degeneracy, i.e. the spindegenerate 
model, where exact results are available as discussed in 
the introduction. As the SPE provide the exact solution 
in the limit of infinite degeneracy one can hardly expect 
reliable results down to Nf — 2. The exact solution for 
Nf — 2 in the symmetric case 2ef + U = and a sym- 
metric half-filled band has particle-hole symmetry, which 
leads to (Po) = (P2) and a symmetric real / - spectral 
function, i.e. the Kondo peak is located at the energy 
zero. In our approach the / and the f 2 subspaces are 
treated in completely different ways and we there cannot 
expect the solution of the SPE for the symmetric Nf = 2 
case to obey particle - hole symmetry. To test how well 
the behavior from the exact solution is reproduced by the 
solution of the SPE we show in fig. |J the probabilities 
(P n ) for Nf = 2, Ef = —2.5 and A = 2.67 as a function of 
U for U > 5. In the exact solution the crossing between 
(Po) and (P2) curves accurs at U = 5, while the SPE 
crossing ist at U — 5.4. If one has symmetric parameters 
2s f + U but Nf > 2 the exact Kondo peak lies at a small 
positive energy and its position for Nf — > 00 is given by 
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the solution of the SPE. While the peak position of the 
exact solution is only defined for integer values of Nt we 
can study the peak position from the SPE as a continu- 
ous function of Nf. This is shown in fig. || for £/ = —2.5, 
C/ = 5 and A = 1.33. With the degeneracy approaching 
Nj = 2 the Kondo peak approaches zero, but does not 
quite reach it. 

The results shown in figs. [| and || indicate that the 
SPE are quite good even for rather small values of the 
degeneracy. A more detailed presentation of results for 
thermodynamic properties and the real / - spectral func- 
tions also at finite temperatures will be presented in a 
forthcoming publication. 

V. SUMMARY 

The generalized mean field theory for the Anderson 
impurity model including double occupancy presented in 
this paper has filled a well known gap in the approximate 
treatments of this model. In contrast to the mean field 
theory for the U = oo model our new approximation can- 
not alternatively be obtained from a simple factorization 
of higher correlations. The use of the coherent state func- 
tional integral and the introduction of two time auxiliary 
fields seems necessary to obtain the correct structure in 
the large degeneracy limit. In the saddle point approxi- 
mation the theory yields the exact result for the ground 
state in the limit Nf — > oo. For large but finite Nf 
the solution of the SPE is expected to provide a good 
approximation for temperatures small compared to the 
Kondo temperature. The spindegenerate limit cannot be 
reproduced exactly, as the theory treats the f° - and/ 2 
- subspaces in completely different ways. The presented 
mean field theory can be generalized in a rather straight- 
forward way from the impurity to the lattice model. 

APPENDIX A: 

In this appendix we present a short but selfcontained 
derivation of the integral equation which has to be solved 
to obtain tie exact ground state in the infinite degen- 
eracy limitEii This equation also describes the exact 
ground state for finite Nf if the valence band is com- 
pletely filled. 

Let | > be the state with all conduction states below 
the Fermi energy, ep = 0, filled and the /-level empty. 
This state couples via Ha to the states 

|£>=4rE^ +c -l 0> ' ( A1 ) 

in which one conduction electron below ef has hopped 
into the /-level. These states couple to states in which 
the /-level is doubly occupied 
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as well as to states with the /-level empty and a par- 
ticle hole-pair in the conduction band. The coupling to 
the latter states is smaller by a factor 1 / y/Nj and can 
therefore be neglected in the large Nf limitEx In this 
limit the ground state takes the form 



b >= A 



> 



/0 />0 f-e 

dea(e) | e > + / de de'b(s,s') | e,e' > 
-B J-B J-B 



where A is the normalization constant. Inserting into 
the Schrodinger equation Ha \ 4>o >= E | O > the 
b(e, e') can be simply expressed-in terms of a(e) and a(e') 
and one obtains the equations!^ 



AE = 



V(e)a(e)de, 



(A4) 



e f - e - AE - (1 - ^-)r(2 £/ +U-e-AE) 



a(s) + V(e) 



Nf 



V{ex)a{e x ) 



B 2e f + U 



AE 



where AE = E - < | H \ >,V(e) = yjN } V{e) 
and 



t{z) 



V{ef 



(A6) 



In the derivation the matrix elements < e \ H \ >= 
V{e) and 

1 



< e,e' | H | e" >= (1 - — ) V(e')S(e - e") + V(s)S(e' - e") 
Nf L 



are used. For the comparison with the saddle point 
equations we define the function c(e) for —B < e < by 

a(e) ee V{e) 



1 + c(e) 



Then equations (A4) and (A5) read 

r° 

AE = 



AE - e f + e + (1 - 1 ^)r(2e / + U-e-AE) 

\V(e) | 2 (l + c(e)) 

AE - s f + e + (1 - ±)f (2e f + U - e - AE) 



c(e) 



J_. 



I V-( £l )| 2 (l + c( £ i)) 



(2E f + U-E-AE-E 1 )AE-£f + Ei + T(2e f + U-ei- AE) 



dsi = ( 



In the limit Ef + U >> B the integral equation ( |A10| ) 
is of separable form. The function c(e) for — B < e < 
can be analytically continued using Eqn. ( |A10| ). Then 
c(z) has a branch cut on the real axis from 2sf + U — AE 
to 2ef + U — AE + B. The saddle point equations reduce 
at T = in the large degeneracy limit to Eqn. (A9) and 
(An]) if one identifies c(e) with x(e) and AE with —A. 
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APPENDIX B: 



In this appendix we derive the expression for the grand 
canonical potential at the saddle point Jsp and discuss 
how it simplifies in the limit Nf — ► oo. 

We begin with the evaluation of the last term on the 
rhs of Eqn. ( |l7| ) for the effective action S e ff. At the 
saddle point the fields F, F, X and X depend on the time 
difference only and the "time-part" of the trace can be 
written as a sum over the fermionic Matsubara frequen- 
cies. If we denote the unperturbed propagators in ac- 
cordance with Eq. (40) as G 00 (iuj n ) we have to evaluate 
tr ln[l — G 00 hi] where the "small trace" runs over the one 
electron Hilbert space. We use tr In a = In det a as the 
determinant is easily calculated using the special form of 

hi © 

det(l - G 00 (iu n )hi(i(j n )) = (1- < / | G 00 | / > VF)(1- <c\G 00 \c>VF) 

- V 2 (b + X) < c | G 00 | c> (6 + X) < f | G 00 | / > 
= l-VF(iu n )G^(iLJ n ) 
(iui n - if)G f f(iuj n ) 



In the second equality we have used (|4lJ) and (|42|) . The 
Matsubara sum is now performed as a contour integral 
using the assumptions about the analytic continuation of 
the auxiliary fields discussed in section II. This yields 



1 





N f Tr In [9 ® l e + l T <8> h + hi] \ 



SP : 



•J A, 



M 



(AJ), 



with 



(AJ)i = N } 



2m 



In 



1 - VF(e + io)G° c °(e + iO) 
l-VF(e-io)G™(e-iO) 



In 



{£-£f-io)Gff{£-io) 
e-e f +io)Gff(e + io) 



and Jjjf + J° the grand canonical potential for a system 
nonintcracting electrons described by Hq. 

If we denote the double integral on the rhs of (|l^) by 
(31 we obtain using the SPE 



I = N f V 



de 
2ttI 



XG 



fc) 



— (FG CC + XGf c ) e +io] 



Using ( p^| ) with Zq(X) ~ 1 for the interesting temper- 
ature regime we finally obtain 



J spa — J A 



M 



(AJ)i -I + X | b | 2 -A (B5) 



For T — > this expression simplifies considerably in 
the large N f limit. Then VF and G ff - G^ are of 



order 1/Nf and the logarithns in (B3) can be expanded. 
To the same order in 1/Nf the Greens functions in (B4) 
can be replaced by the unperturbed functions Gj?° and 
G°j c . Then in order (1/Nf)° all terms in the expression 
for J spa cancel except the first and the last one. With 
/4 0) =< | H | > this yields 



(E 



,SPA 



= E> 



(0) 



(B6) 



as mentioned following Eqn. 
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APPENDIX C: 



In this appendix we discuss how the "real /- 
occupancy" in the saddle point approximation {nf)^ a f 
can be expressed in terms of the "pseudo-/-occupancy" 
( n f)ps A = N f GS ff A ( T = -°)- The calculation shows 
that the constant A defined following Eq. (|34|) presents 
the probability for the double occupancy of the f -level. 

The real /-occupancy can be calculated from the grand 
canonical potential J = — In Z/(3 with Z defined in (||) 
by differentiation with respect to the /-level energy ej 

dJ 

(rif)real = (CI) 
ClBf 

In the SPA the grand canonical potential is given in 



(|B5|). With h = ho + hi it reads 

Jspa = —(Nf/f3)Trln(d + h)-I + \(\b\ 2 -1) 

= J h -I + X(\ b | 2 -1) (C2) 

In the following the label SPA at all fields and 
Greens functions will be suppressed. The contribution 
to (rif) rea i from the first term on the rhs of (C.2) is 

57 " ev*^f;> (cs) 



and with h = l r <S> ho + hi given in ( |18| - |19| ) we obtain 

^ = (1 + ^){n f)ps + - [Nf Vtr T{ G ff -) + N f Vtr r{ G cc -) 
+ NjV^-GfM + §-G cf (0)) + N f Vtr T (G fc ^- + G c/ ||)] 

where tr T denotes the "time" -integrati on p art of the 
trace. Differentiating the double integral (B4) yields 

— = - / drdr'-^- L >_ x{t _ t i )+x[t _ t i ) y 

d£f PJ G u h {T-T') L de f de f 

+ ^ (T ., )+ , (T _, ) ^ ] _ (2+ - ) , 

and from the last term on the rhs of (C.2) we obtain 

^ m .. 1 ,_« (l „.. 1 , +A( « t+ ,» ) 

(C6) 



If we add up all contributions and use the SPE ( J27| - |32| ) 
and ( p^ ) we obtain the simple result 

(n f ) rea i = (n f ) ps + 2A (C7) 

This suggests that A is the probablity for the double 
occupancy of the /-level. This is readily verified by cal- 
culating < P2 >, where P n is projection operator for 
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the /-level occupied by n electrons. A calculation very 
similar to the one presented above yields in SPA 

<P2>=^j=A (08) 

In the limit Nf — > oo and T — > this identity can be 
explicitely checked using the results of appendix A and 
the expression for < P-2 > presented in ref. 13. From 
(C.7) and (C.8) and J2n=o < -^n >= 1 we also obtain 
< Pi >= (nf)ps and < Pq >=| b | 2 , which elucidates the 
meaning of the SPE (Q . 

The generalization of (C.8) beyond the SPA is given 
by using ( |l6|) to evaluate dJ/dU. This yields 



where we have taken Zf}(X) ~ 1 as usual. If we take 
the FI on the rhs of (C.9) as the generalized expression 
for A, equation (C.7) also holds generally. 



1 P.W. Anderson, Phys. Rev. 124, 41 (1961) 

2 P.A. Lee, T.M. Rice, J.W. Serene, L.J. Sham, and J.W. 
Wilkins, Comments on Condensed Matter Physics XII, 
No. 3 (1986) 

3 O. Gunnarsson and K. Schonhammer, in: Handbook on 
the Physics and Chemistry of Rare-Earths, K. Gschnei- 
dner, Jr., L. Eyring, and S. Hiifner eds. (North-Holland, 
Amsterdam, 1987), Vol. 10, p. 103 

4 N.E. Bickers, Rev. of Modern Physics 59, 845 (1987) 

5 P. Schlottmann, Physics Reports 181, 1 (1989) 

6 O. Gunnarsson and K. Schonhammer, Phys. Rev. Lett. 50, 
604 (1983) 

7 O. Gunnarsson and K. Schonhammer, Phys. Rev. B28, 
4315 (1983) 

8 Y. Kuramoto, Z. Phys. B53, 37 (1983) 

9 P. Coleman, Phys. Rev. B29, 3035 (1984) 

10 N. Read and D.M. Newns, J. Phys. C16, 3273, ibid. L1055 
(1983) 

11 A. Houghton, N. Read, and H. Won, Phys. Rev. B35, 5123 
(1987) 

12 P. Coleman, Phys. Rev. B35, 5072 (1987) 

O. Gunnarsson and K. Schonhammer, Phys. Rev. B31, 
4815 (1985) 

14 J. Holm and K. Schonhammer, Solid State Commun. 69, 
969 (1989) 

15 Th. Pruschke and N. Grewe, Z. Phys. B74, 439 (1989) 

16 M. Jarrel, Phys. Rev. Lett. 69, 168 (1992) 

A. George and G. Kotliar, Phys. Rev. B45, 6479 (1992) 

17 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989) 

18 M. Jarrell, Phys. Rev. Lett. 69, 168 (1992) 

19 Th. Pruschke, D.L. Cox, and M. Jarrell, to be published 

20 CM. Varma and Y. Yafet, Phys. Rev. B13, 2950 (1976) 



18 



Z. Zou and P.W. Anderson, Phys. Rev. B37, 627 (1988) 
J.H. Kim, K. Levin, and A. Auerbach, Phys. Rev. B39, 
11633 (1989) 

H. Kaga and T. Saikawa, Physica C172, 400 (1991) 



FIG. 1. The spectral functions pp, p F and px for Nf — 14, 
e f = -2.5, U = 8 and A = 2.67. 

FIG. 2. The pseudo /-spectral function {pff) ps - The in- 
set shows the the additional spectral weight in the interval 
[iu,£u + B\. The parameters are the same as in fig. 1 

FIG. 3. pcc(e) and p/ c (e) are shown. The parameters are 
the same as in fig. 1 

The spectral functions 

FIG. 4. The occupation probabilities (P n ) for Nf — 2, 
Sf — —2.5 and A = 2.67 as a function of U. 



FIG. 5. The position of the Kondo peak as a function of 
Nf for s f = -2.5, (7 = 5 and A = 1.33. 
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